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Abstract. In this paper we investigate the energy distribution of states of 
a linear scalar quantum field with arbitrary curvature coupling on a curved 
C*^ ' spacetime which fulfill some local thermality condition. We find that this con- 

^ i dition implies a quantum energy inequality for these states, where the (lower) 

, energy bounds depend only on the local temperature distribution and are lo- 

cal and covariant (the dependence of the bounds other than on temperature 
in - is on parameters defining the quantum field model, and on local quantities 

constructed from the spacetime metric). Moreover, we also establish the av- 
eraged null energy condition (ANEC) for such locally thermal states, under 
growth conditions on their local temperature and under conditions on the free 
parameters entering the definition of the renormalized stress-energy tensor. 
These results hold for a range of curvature couplings including the cases of 
$"H , conformally coupled and minimally coupled scalar field. 
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1. Introduction 



One of the problems in quantum field theory is that the space of states is enor- 
mous, and that it is difficult to establish criteria which single out the states of in- 
terest in various physical situations. While this problem is met only in a mild form 
when considering quantum fields on Minkowski spacetime, it becomes quite pressing 
' when trying to combine quantum field theory and gravity. A situation where this is 

predominant is quantum field theory in curved spacetime, where quantized matter 
fields propagate on a classical background spacetime which may be curved, and 
where the spacetime metric is not stationary - spacetimes of Friedman- Robertson- 
^ ' Walker type, for instance, are of particular interest in this context as they are simple 

models for cosmological scenarios. In trying to model the conditions of the early 
stages of the universe, one would like to distinguish quantum field states which 
are, at least locally, not too far from thermal equilibrium, and for which one can 
assign, at least locally, a temperature. It is not at all easy to arrive at a meaningful 
concept of temperature for states of a quantum field on a generic spacetime since 
the (global) notion of temperature makes reference to a global inertial frame, and 
that is not available in the presence of spacetime curvature. Another manifesta- 
tion of this circumstance is the observer-dependence of the concepts of particle and 
temperature, as illustrated by the Fulling-Unruh, and related effects [571 12"51 12"4"] . 

Nevertheless, following a proposal by Buchholz, Ojima and Roos [5] and further 
investigated in [3 [TO], it is possible to introduce a covariant concept of states which, 
at given points x in spacetime, look like thermal equilibrium states with respect 
to a certain set, S x , of reference-observables. Such states are called local thermal 
equilibrium (LTE) states with respect to S x . While there is some leeway in the 
determination of S x , it is important that S x does not contain observables which are 
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sensitive to flux-like quantities, but instead observables which correspond, in the 
situation of global thermal equilibrium states defined with respect to an inertial 
frame (in the absence of spacetime curvature) , to intensive thermal quantities. For 
a (scalar) quantum field <fi(x), typical elements of S x are the Wick-square : 4> 2 : (x) 
of the field and its so-called "balanced derivatives" . We will discuss the concept of 
LTE states further in the next section. However, it is important to mention that for 
an LTE state u> of the massless linear scalar field on a generic (globally hyperbolic) 
spacetime, the expectation value of the Wick-square of the field at any spacetime 
point x, 

(1) T 2 (x) = (: (f) 2 : (x))^ = lo{: 4> 2 : (x)) 

equals, up to a constant, the square of the absolute temperature of the state at x. 
(A similar statement holds for the linear scalar field with positive mass parameter.) 

Quite clearly, the local thermodynamic properties of states are linked with the 
local energetic properties of states, and investigation of that relation is the topic 
of the present article. A quantity of prominent interest in quantum field theory in 
curved spacetime - and with relevance to questions in cosmology - is the expecta- 
tion value of the stress-energy tensor, (T a f,(a;)) w , in a state uj for a quantum field 
on a curved spacetime. It becomes particularly important when considering the 
semiclassical Einstein equations (in geometric units), 

(2) G ab {x) = 87r(T a t ss (x) + (T ab (x)) u ) , 

where G ab (x) is the Einstein tensor of the spacetime geometry at spacetime point x 
and T£ l b ass is the stress-energy tensor of macroscopically modelled (classical) mat- 
ter. If the gravitational curvature effects caused by the macroscopic energy/matter 
distribution are very high, this may induce quantum field theoretical "particle cre- 
ation effects" , as a result of which (T a ;,(a;)) w may turn out as a significant correction 
to the macroscopic matter distribution, depending on the quantum state to. One 
may expect that this is a realistic scenario as nowadays the Casimir force puts limits 
to the design of devices in micro- and nano-technology [351 II] > an d this effect can 
be seen in a similar vein. 

One of the interesting features of the expectation value of stress energy is that the 
energy density seen by an observer travelling on a timelike geodesic 7 with tangent 
vector v a at x, (T ab (x)} LJ v a v b , is unbounded above and below as lo ranges over the 
set of all states ui (for which the expectation value of stress-energy at any spacetime 
point x can be reasonably defined) . This is a long known feature of quantum field 
theory (see and is in contrast to the behaviour of macroscopic matter which 
can - with good motivation - usually be assumed to satisfy one of the classical 
energy conditions, like the weak energy condition, which means T^ l b ass (x)v a v b > 0, 
i.e. the energy density seen by any observer is always positive at any spacetime 
point x. 

Energy positivity conditions like the (pointwise) weak energy conditions play an 
important role in the derivation of singularity theorems [251 40]. One consequence 
of energy positivity conditions when plugged into Einstein's equations is that grav- 
itational interaction is always attractive. Negative energies, in contrast, would be 
affected by a repelling gravitational interaction. This could, a priori, lead to solu- 
tions of Einstein's equations exhibiting very strange spacetime geometries, such as 
spacetimes with closed timelike curves, wormholes or "warpdrive scenarios" |30l fT]. 
Moreover, concentration of a vast amount of negative energies and their persistence 
over a long duration could lead to violations of the second law of thermodynamics. 

Motivated by the latter point, L. Ford has proposed that physical states of quan- 
tum fields in generic spacetimes should not permit arbitrary concentration of large 
amounts of negative energy over a long duration [5T] . Such limitations on physical 
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quantum field states have come to be called quantum energy inequalities (QEIs). 
Let us explain this concept in greater detail. Suppose that 4>{x) is a quantum 
field on a generic spacetime. (Actually, only smeared quantum field quantities like 
4>(f) = J f(x)4>(x) dvo\(x) define proper quantum field operators operators; dvol(x) 
denotes the metric induced spacetime volume form and / a smooth, compactly sup- 
ported test-function. The quantum field can be of general spinor- or tensor type, 
but we suppress any corresponding indices here.) Then let £ be a set of states 
of the quantum field such that the expectation value of the stress-energy tensor, 
(T a h(a;)) w , is defined for each w £ § at each spacetime point x. We will furthermore 
suppose that this quantity is continuous in x for each lu G £. Under these assump- 
tions, we say that the set of quantum field states L fulfills a QEI with respect to 7 



holds for all smooth (or at least C 2 ) real functions h having compact support on 
the (open) curve domain, with a constant 9(7, h) > — 00; the constant may depend 
on the curve 7 and the weighting function h, but is required to be independent of 
the choice of state uj e £. 

There is a limiting case of a QEI: If 7 is a complete (lightlike or null) geodesic, 
then one says that a set of states L fulfills the averaged null energy condition 



holds for all states ui £ L. Conditions of such form (and related conditions, see 
if valid for all complete null geodesies, allow conclusions about focussing 
of null geodesies for solutions to the semiclassical Einstein equations similar to 
that resulting from a pointwise null energy condition [36] [32] [42] . (See also the 
beginning of Section 4.) Thus, the ANEC is a key property for deriving singularity 
theorems for solutions to the semiclassical Einstein equations. 

Quantum energy inequalities have been investigated extensively for quantum 
fields subject to linear field equations in the recent years, and there is now a wealth 
of results in this regard. We refer to the reviews by Fewster and by Roman [33] fT3j 
for representative lists of references. Important to mention, however, is the fact that 
for many linear fields, like the minimally coupled scalar field, the Dirac field and 
the electromagnetic field, it could be shown that the set of Hadamard states fulfills 
a QEI with respect to timelike curves 7 in generic globally hyperbolic spacetimes 
P~2l QjH [15] . Hadamard states are regarded as physical states in quantum field the- 
ory in curved spacetime, and expectation values of the stress-energy tensor at any 
given spacetime point are well-defined for these states (up to finite renormalization 
ambiguities), cf. [41] for discussion. There is also an intimate relation between 
QEIs, the Hadamard condition and thermodynamic properties of linear quantum 
fields [20 . It has been shown that QEIs put strong limitations on the possibility 
of solutions to the semiclassical Einstein's equations to allow exotic spacetime sce- 
narios such as wormholes or warpdrive [221 1171 131] . It is also worth mentioning 
two other recent results. First, it has been shown that the non-minimally coupled 
linear scalar field on any spacetime violates QEIs for the class of Hadamard states; 
nevertheless, the class of Hadamard states fulfills in this case weaker bounds, called 
"relative QEIs", cf. [14] for results and discussion. Secondly, one is interested in 
lower bounds q-y(h) which depend (apart from renormalization constants entering 
the definition of expectation value of the stress energy tensor) only on the under- 
lying spacetime geometry in a local and covariant manner, and one also aims at 
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making this dependence as explicit as possible. Considerable progress on this is- 
sue, for the case of the minimally coupled linear scalar field on globally hyperbolic 
spacetimes, has been achieved in [18] . 

In the present article, we will derive QEI-like bounds on sets of LTE-states of 
the non-minimally coupled linear scalar field <j>{x) on generic globally hyperbolic 
spacetimes. More precisely, we consider LTE states u> whose thermal function 
^(x) = (: 4> 2 : is bounded by some constant Tq (corresponding to a maximal 

squared temperature) and we will show that there are upper and lower lower bounds 
for the averaged energy density 



averaged against a C 2 -weighting function r\ > with compact support along any 
causal geodesic 7 with afhne parameter r and tangent v a = 7°. The lower bound 
depends only on Tq, the geodesic 7 and rj, while the upper bound depends addition- 
ally on local tetrads entering into the definition of LTE states. The lower bound is 
therefore state-independent within each set of LTE states u> with a fixed maximal 
value of $ w . The bounds depend on the spacetime geometry in a local covariant 
manner which, together with their dependence on Tq, we will make explicit. This 
result holds for all values of curvature coupling £ in the field equation @, and 
upon averaging along causal geodesic, not only those which are timelike. Hence, 
the result is not immediate from know quantum energy inequalities for Hadamard 
states, as these are violated in general for non-minimally coupled fields [14j . and 
upon averaging along null geodesies [15] . Furthermore, we will show that the ANEC 
holds for LTE states to of the quantized linear scalar field with curvature couplings 
< £ < 1/4, provided that the growth of the thermal function $ w along the null 
geodesies 7 fulfills certain bounds. Despite the fact that we have to assume that 
the LTE states we consider are Hadamard states - in order to have a well-defined, 
local covariant expression of expected stress-energy for these states - our derivation 
of QEIs and ANEC makes no further use of the Hadamard property but uses only 
properties of LTE states. Therefore, one may expect that, in principle, similar re- 
sults could be derived for LTE states of interacting quantum fields. This prospect 
can actually be seen as one of our motivations in view of the fact that quantum en- 
ergy inequalities seem to be very difficult to obtain (if valid at all) for very general 
sets of states in interacting quantum field theory, and that, on the other hand, one 
may argue that only special classes of states are of physical interest. We will come 
back to this point in Sec. 6. 

This article is organized as follows. We will discuss the concept of LTE states, as 
far as needed for our purposes, in Section 2. In Section 3 we derive upper and lower 
bounds for the geodesically averaged expectation values of energy density for LTE 
states. The validity of ANEC for certain LTE states will be studied in Section 4. 
In Sec. 5 we indicate that the results of Sees. 3 and 4 hold also for a more general 
notion of LTE states. We conclude with discussion and outlook in Sec. 6. 



The system under investigation in the present article is the non-minimally cou- 
pled linear scalar field on globally hyperbolic spacetimes. A globally hyperbolic 
spacetime will be denoted by a pair (M, g) where M is the spacetime manifold 
(assumed to be C°°) and g is the Lorentzian metric. We will consider the case 

of spacetime dimension equal to 4 with metric signature (H ), but most of 

our considerations can be readily generalized, with appropriate modifications, to 
arbitrary spacetime dimensions. We recall that global hypcrbolicity means that the 




2. Local Thermal Equilibrium States 



spacetime is time-orientable and possesses Cauchy surfaces [40 , 2] . Our conventions 
for curvature quantities, like in [14 , are those of Birrell and Davies, i.e. [-,-,-] in 
the classification scheme of Misner, Thorne and Wheeler. 

The classical linear scalar field jjona globally hyperbolic spacetime (M, g) obeys 
the field equation 

(4) (V^V M +£i? + m 2 V = 

where V is the covariant derivative of g and R is the scalar curvature corresponding 
to <?; the constants £ > and m > are the curvature coupling and the mass 
parameters, respectively. The case £ = corresponds to minimal coupling. 

The quantization of the system proceeds as follows. Owing to global hyper- 
bolicity, there are (for each fixed £ and m) two uniquely determined linear maps 
E ± : Cg°(M,R) -> C°°(M,R) so that 

(5) E ± (WV tl + £R + m 2 )f = f= (V^V M + £R + m 2 )E ± f 

holds for all / G Cq° (M, R) , and additionally, supp( J B ± /) c J ± (supp(/)), where 
J ± (G) is the causal future/past set of G c M [2 J. These are called the ad- 
vanced/retarded fundamental solutions of the wave-operator (V^V^ + £i? + m 2 ), 
and with their help one can construct the real bilinear form 

(6) /a) - / fi(x)(E~f 2 - E+f 2 )(x) dvol(x) 

J M 

on C§°(M, E) which turns out to be antisymmetric. Note that S is uniquely de- 
termined by (M,g), £ and m. Fixing £ and m, one can now define the complex 
^-algebra A(M,g) = A((M,g),£, m) with unit element 1 as being generated by a 
family of objects <f>(f), f G C£°(M,R) which are required to fulfil the following 
relations: 

(a) / i—* </>(/) is real-linear, 

(b) m* = <m/), 

(c) 0((VV M +£i? + m 2 )/) =0, 

(d) [0(A), 0(/ 2 )] = i#(/i,/ 2 )l. 

Here [A, B] = AS — denotes the commutator. Since the generators 4>(f) of 
7l(M, 5) obey, according to (d), the canonical commutation relations in a covariant 
manner, one has thus obtained a quantization of the system in an abstract form. 
The hermitean elements in A{M, g) correspond to observables of the quantized sys- 
tem, but they do not contain all observables that one may wish to consider, so that 
the algebra A(M, g) will have to be enlarged to include those additional observables 
as well. We will come back to this point. For the moment, a state u> of the quan- 
tized linear scalar field is, by definition, a linear functional lo : A{M,g) — > C, 
A 1 ► = (A) u , with the additional property that u> is positive, meaning 

io(A*A) > for all A G A(M,g), and also with the property that u> is normal- 
ized, i.e. w(l) = 1. Now it is known from examples that not every state according 
to this definition corresponds to a physically reasonable configuration of the system 
and that selection criteria for physical states are needed. In the case of the linear 
fields on curved spacetime, the best candidates for physical states are quasifree 
Hadamard states, and most other physical states can be derived from those 38J. 

We will very briefly summarize the concept of a quasifree Hadamard state. (For 
a more in-depth discussion, see [IE]-) For any state u> on A, the n-point functions 
are the maps 



(7) 



fn) = U,M/i) • • • <K/n)) • 
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Clearly, each state on A(M,g) is determined by all the n-point functions. A 
quasifree state to on A(M,g) is a state which is entirely determined by its two- 
point function, by requiring that the truncated n-point functions vanish [6]. This 
can also be expressed as oj{e lt ^^) — e~* *2"(/>/)/ 2 ; to be interpreted as a se- 
quence of relations in the sense of formal power series in t. A (quasifree) state lo is 
called Hadamard state if its two-point function is of Hadamard form. This is the 
case, in turn, if for any geodesic convex neighbourhood N of any given point x , 
and any time function t on the underlying spacetime M, one can find a sequence 
H% eC k (N x N, C) such that for all /i,/ 2 € C* °°(A,R) one has 

(8) *Wi,/ 2 ) 



lim -L / (G kie {x,x')+H k (x,x , ))f 1 (x)f 2 (x , )dvol(x)dvol{x'), 
' *0+ Air 2 J NxN 



where 

(9) Gk, e (x,x') = U{X ' X> 



a{x, x') + 2i(t(x) - t(x'))e + e 2 
+V fe (x, x')\n(a{x, x') + 2i(t(x) - t{x'))e + e 2 ) . 



Here. 



(10) V k (x,x') = ^U j (x > x')a(x,x'y , 

3=0 

a(x,x') is the squared geodesic distance from x to xQ, and U and Uj are smooth 
functions on TV x N determined by the Hadamard recursion relations. Thus, the 
term G k ,e is, for each fc, determined by the local spacetime geometry and the 
parameters £ and m of the scalar field equation. For later use, we define the 
distribution 

(11) Wi,h) = £ lim ^2 / G k>£ (x,x')f 1 (x)f 2 (x')dvol(x)dvol(x') 

for test-functions f\, f% supported in a geodesic convex neighbourhood N. 

It is worth noting that the existence of very many quasifree Hadamard states 
(spanning an infinite dimensional space) has been established for the linear scalar 
field on all globally hyperbolic spacetimes. Moreover, in stationary, globally hy- 
perbolic spacetimes, the canonical ground state as well as the thermal equilibrium 
states (KMS states) are known to be quasifree Hadamard states [34j . 

Let us now turn to the concept of local thermal equilibrium states introduced in 
[S] and further investigated in [51 [TD] . This will be done first for the case that the 
underlying spacetime is just Minkowski spacetime. Our discussion here is limited 
to the linear scalar field, but as explained in [jj], the discussion can be general- 
ized to include general quantum field theories. For the quantized linear scalar field 
on Minkowski spacetime, there is in each Lorentz frame a unique quasifree ther- 
mal equilibrium state at given temperature. Actually, for fixed temperature this 
state depends only on the time-direction of the Lorentz frame. Let eo be such 
a time-direction, i.e. a timelike, future-pointing unit vector on Minkowski space- 
time, and let e%, e?, e% be a set of spacelike unit vectors so that e = (eo, e%, e%, e^) 
forms an orthonormal tetrad on Minkowski spacetime. When choosing coordinates 
(x , x , x 2 , x 3 ) on Minkowski spacetime such that the coordinate axes are aligned 
with the tetrad, the two-point function of the unique quasifree thermal equi- 
librium (KMS) state w /3e at inverse temperature /3 > 00 with respect to the Lorentz 



1 Following [IBj, we choose a positive for x and x' spacelike related and negative for x and x' 
timclike related, so that e.g. on Minkowski spacetime it is given by <j(x, x') = —g a i,(x—x') a (x—x') b 
2 f3 = l/(kgT) where kg is Boltzmann's constant and T is absolute temperature 
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frame defined by e is given by 

d 4 p 



(12) Wt{x, x') = / e-^-^^e{ Po )5{p^ - m 2 ) 



(2tt) 3 (1 - e-/ 3 ^) ' 

to be interpreted in the sense of distributions, where e(po) is the sign function of 
Pq. The uniqueness implies that every intensive thermal property (e.g. pressure, 
density etc) can be expressed as a function of the timelike vector j3eo. 

The passage from these global equilibrium states to states of local thermody- 
namic equilibrium now uses spaces S x of observables located at spacetime points 
x. Mathematically this implies that these observables are no longer defined as 
operators but only as quadratic forms (their products are usually not defined). In 
physical terms, the observables in S x should model idealized limits of measurements 
of intensive thermal properties of states in smaller and smaller spacetime regions. 
This culminates in the requirement that, for s{x) £ S x , the functions 

have to be independent of x and non-constant as functions of j3. Furthermore, for 
many s(x) G S x one can identify the thermal quantity to which $ s , as a function 
of (3, actually corresponds. As an example, one calculates that for s(x) =: <fi 2 : (x) 
(Wick-square) one obtains for the massless case 

1 k 2 T 2 

which leads to the identification of this observable as a "scalar thermometer" , giv- 
ing the square of the local temperature T times some fixed constant. Choosing 
appropriate units to measure T, this constant can be set equal to 1, justifying the 
notation T 2 for the expectation value of the Wick-squarc in thermal equilibrium 
states already alluded to above. Actually this seems quite similar to what one 
would do to construct a thermometer in the laboratory: Take some (small) device 
which, when exposed to a situation known to be in equilibrium at some temperature 
T, gives a reading which is a simple function of T. 

In the investigations of local thermal equilibrium for linear scalar fields 4>{x) 
on Minkowski spacetime carried out in previous articles HH] , the spaces S x 
are chosen as the linear spaces generated by elements s(x) = 3^ : 4> 2 : (x), re- 
ferred to as the balanced derivatives of the Wick-squared field : <fi 2 : (x). |j Here, 
fi = (/ii, . . . , /i n ) € N™ is a multi-index of arbitrary length n, and the balanced 
derivatives are defined by 

(13) &V : 4? : (x) = Hm 0„ {4>{x + ()0(x - C) - ^(^{x + Qcj>(x - 0)1) 

where <9 M = d^i ■ ■ ■ d^>i„ , and where uj vac is the vacuum state. The limit is taken 
along spacelike directions C, so that <f>(x + ()(f>{x — Q is well defined as a quadratic 
form, and the limit defines an operator- valued distribution after smearing in x with 
test-functions. For multi-index length equal to 0, the balanced derivative equals just 
the Wick-square : <f> 2 : (x). For linear fields on Minkowski spacetime, this definition 
of the Wick-square coincides with the usual normal ordering prescription. Owing 
to the translation invariance of the KMS-states uj^ e one can easily check that the 
thermal functions $9^2. = w' 3e (3 M : (f> 2 : (x)) are independent of x. 

Following [H [S] , a state of the linear scalar field on Minkowski spacetime is said 
to be locally in thermal equilibrium at a spacetime point x if it looks like a global 
thermal equilibrium state uj@ e as far as the expectation values of elements in S x are 
concerned. The following definition, taken from [9], expresses this more formally. 



^In the references [9]|8], the notation 9** is used, but we prefer to view fj, as a co-tensor index 
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2.1. Definition. A state uj of the quantized linear scalar field <ft(x) on Minkowski 
spacetime is called S^-thermal at the spacetime point x if there are an orthonormal 
tetrad e with eo timelike and future-pointing, and (3 > 0, such that 

(14) lj(s(x)) = uj^ e (s(x)) 

holds for all s(x) £ S x , where S x is spanned by 3^ : (ft 2 : [x) as /i ranges over all 
multi-indices. 

If is some open set of spacetime points, a state to of the quantized linear scalar 
field on Minkowski spacetime is called So-thermal if (|14[) holds and u>(s(x)) varies 
continuously with x, for s(x) £ S x and x € 0. That means, lu is S x -thermal at 
each x £ 0, where (3 and e in (| 14[) may vary with x. 

For So-thermal states, the expectation value (^ s ) UJ (x) of an extensive thermal 
quantity <& s at x, whose local measurement is modelled by s(x) £ S x , is then 
given as ($ s ) w (x) = v(s(x)). This leads for So-thermal states to an assignment 
of thermal quantities to each x S whose values in general vary with x, and this 
assignment is consistent in the sense that relations among the thermal quantities 
(like equations of state) also hold at each point. In the case of s(x) —: (ft 2 : (x), 
one obtains in this way for an Sq -thermal state oj an assignment of the expected 
squared temperature to each spacetime point x £ O. 

The concept of states which are locally in thermal equilibrium has also been 
generalized in [9] to allow mixtures of global thermal equilibrium states on the 
right hand side of (fT4|) . We will summarize this generalized concept in Sec. 5. 

In attempting to extend the concept of local thermal equilibrium states to quan- 
tum fields in curved spacetime, one faces a couple of difficulties which are, of course, 
connected to the occurrence of curvature and the related lack of global vacuum 
states and global equilibrium states. Primarily, these difficulties are: 

(i) The definition (fl3|) of balanced derivatives 3 M : (ft 2 : (x) uses the Minkowski 
vacuum state w vac as preferred vacuum state. 

(ii) Moreover, the definition (fT5|) uses the affine space structure of Minkowski space- 
time. 

(iii) Def. 12.11 uses global thermal equilibrium states ui^ e on Minkowski spacetime, 
for which there is no counterpart on generic curved spacetimes. 

Thus, there is no verbatim translation of the concept of local thermal equilibrium 
states given in Def. 12.11 

It is clear that problems (i) and (ii) concern the definition of balanced derivatives 
of a Wick-squared quantum field in curved spacetime. We will soon turn to that 
problem. Assuming that the definition of balanced derivatives in curved spacetime 
is settled, a proposal was made in [10] to surpass problem (iii) . The idea is to define 
that a state u> of the quantized linear scalar field (ft on a curved spacetime (M, g) is 
S x -thermal at a point x in M iiw(s(x)) = ujP e (s (x )) holds for all s(x) £ S x . Here, 
Lu@ e is a thermal equilibrium state of the free scalar field <f> (with same parameters 
as (ft) on Minkowski spacetime M a , x is a point in M , and s (x ) is the flat 
space counterpart of s(x). To explain what this latter phrase means precisely can 
be seen as part of the definition of balanced derivatives in curved spacetime, but 
certainly one would require that s (x ) corresponds to a balanced derivative of the 
Wick-square of <f> Q if s(x) corresponds to a balanced derivative of the Wick-square 
of (ft. In the approach of [10] . the requirement of local thermality on u> is thus not 
implemented by comparing expectations values of pointlike thermal observables 
with the corresponding expectation values in a global thermal state (as such states 
need not exist) , but with the "flat space version" of a thermal equilibrium situation 
for the quantum field. The motivation for this approach is that S x -thermality is a 
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pointwisc property which should not be affected by curvature; this, in turn, rests 
largely on the equivalence principle. 

For a linear scalar quantum field </> on a curved spacetime it is simple enough 
to know what its flat space counterpart <j) should be. However, as pointed out in 
|10j , one may invoke the concept of a local covariant quantum field theory [JJ 126) 
to know this also for more general types of quantum fields. The concept of local 
covariance affects also the elements s(x) G S x ; in our situation where we start from 
a linear scalar field <f> in curved spacetime - which is known to have the structure 
of a local covariant quantum field - the balanced derivatives of Wick-squares of <f> 
should be defined in such a way that they are also local covariant quantum fields. 

Let us thus discuss our proposal for the generalization of the concept of balanced 
derivatives in a curved spacetime. Our discussion is greatly facilitated by the cir- 
cumstance that for the purpose of deriving QEIs for LTE states we need only focus 
on balanced derivatives up to second order, corresponding to a multi-index length 
of (J, not greater than two. Accordingly, we will define the LTE property on curved 
spacetime only with balanced derivatives of the Wick-squared field up to second 
order, see below. 

We proceed in two steps. First, we shall consider the generalization of expressions 
like d^nd^v f(%+C, x ~ C)lc=o for C 2 -functions / from Minkowski spacetime to curved 
spacetime. This discussion is entirely of differential geometric nature. In a second 
step, we have to give a generalization of the quantity 

f(x + c,x-Q = cj(cf>(x + $4>{x - 0) - ^ vac W>(x + CM* - 0) 

for Hadamard states u> on curved spacetime where there is no counterpart of w vac . In 
doing this we have to ensure, as mentioned, that the resulting balanced derivatives 
of the Wick-ordered linear scalar field give rise to local covariant quantum fields. 

Turning to the first step, let (M,g) be a spacetime and suppose that N is a 
geodesically convex neighbourhood of some point x in M. The exponential map 
at x will be denoted by exp^,. A fairly obvious generalization of the first balanced 
derivative of a function / S C 2 (N x N) arises by requiring 

f(exp x (Xv),exp x (-Xv)) 

for all spacelike vectors v = v a e T X M lying in exp~ 1 (iV). By linearity, this 
determines a co- vector 3 Q /(x) in T*M. We define T{ r s r s ',)(Al x M) as the bundle 
over M x M whose fibre at {y,y') £ M x M is given by T y ( r s )M® T y ,( r s ',)M, where 
T y { r s )M coincides with the space of r-fold contravariant and s-fold covariant tensors 
at y. If (y,y r ) i— » V(y,y') is any C 1 section in T{ r a r a ,)(M x M), then we denote by 
V ' oV the covariant derivative with respect to the y-entry and by V a >V the covariant 
derivative with respect to the y'-entry. Furthermore, we denote by I^Lz— V(x,x) 
the coincidence value of V(y,y') for y — x = y' . With these conventions, one has 

9„/(a:) = V„/U-V.'/L», 

and if / is C 2 , this defines a C 1 co- vector field as x varies. 

The second order balanced derivative B a bf(x) can then be defined as follows. 
V a f(y,y') is a y'-dependent co-vector at y, and vice versa for V a 'f(y,y'). Let 
v = v a be a (spacelike) vector in T X M which lies in exp~ 1 (N), so that rj v : A i— * 
exp x (Xv) (—1 < A < 1) is the geodesic determined by v at x. Correspondingly, we 
can define the map of parallel transport P V: \ : T c * xp <\ V \M — > T*M of co- vectors 
from exp A „ = T) V (X) to x = rj v (0) along the geodesic rj v . A geometrically natural 
definition of the second order balanced derivative <5 a bf(x) of / at x is then obtained 



V a 3af(x) = 4t 
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by demanding that 

d 



v a w b B ab f(x) = 

d\ 



d 



dX 



w P_„^V b '/(exp :z ,(Aw),exp a ,(-Aw)) 

A=0 

holds for all (spacelike) vectors v,w g T X M with v G exp~ 1 (7V). Using the proper- 
ties of the parallel transport, it follows that M 3 x <3 a bf{x) is a continuous (if / 
is C 2 ) (^-tensor field on M, and 

(15) d ab f(x) = V a V fe /U-V W/U 

-v 'V 6 /Lx+v a 'V 6 '/Lx- 

Note here that the covariant derivatives on the right hand side act on y for unprimed 
indices and on y' for primed indices, and primed and unprimed tensor indices are 
identified at the coincidence point y — x = x' . 

Turning to the second step, suppose that (M, g) is a globally hyperbolic space- 
time, and that N is a geodesic convex neighbourhood of some point x E M . Then 
define the distributions (k > 2) 

Wi,/a) = /a) + mh, h)) + i*(fi,h) , h,h G C?(N,R), 

where % and $ are defined above for the quantized linear scalar field on (M, g). 

Next, let w be a quasifree Hadamard state of <j> on (M, g); then define the point- 
split renormalized two-point function obtained by subtracting the symmetrized 
Hadamard parametrix (SHP) % from the two-point function: 

^T(/i, h) = n'ifu h) - %(h, /a) , /i,/ 2 e C °°(iV, R) . 

A first observation is that Wjf is symmetric, #^J p (/i,/ 2 ) = ^, S | P (/a,/a)- Fur- 
thermore, for k > 2, ^ S fc P is given 

^T(/i, /a) = / ^ H fc p (y, y')fi(y)f 2 (y') dvo\(y) dvolfo') , /i,/ 2 e C °°(7V, K) , 

JjVxA' 

with W^ 1 ^ G C 2 (N x AT) for k > 2, which follows from the definition of Hadamard 
form (see also the arguments in [IB])- 

Consequently, one may now define for any quasifree Hadamard state u> the expec- 
tation value of the SHP- Wick square of (j> and the corresponding second balanced 
derivatives in the following way: 

2.2. Definition. 

(16) w(: 2 :shp (x)) = W^[ x , 

(17) *:^: S hpW) = dabW^flcc 
for x S M, lOT'tt fc > 2. 

We add a few observations to this definition. 

(a) One can likewise define the first balanced derivative of u>(: 4> 2 :shp (%)), but 
since W^ p is symmetric, its first balanced derivative vanishes. This is similar 
to the property of balanced derivatives of the Wick-square of the quantized linear 
scalar field <p on Minkowski spacetime, which can be traced back to the symmetry 
of 

u; (My)My'))-u vac (My)My')) 

with respect to y and y' , for each quasifree Hadamard state uj of <p . This provides 
motivation why we define the Wick-ordering by subtraction of the symmetrized 
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Hadamard parametrix. 

(/3) Actually, W^F^y, y') and V a V{,' W^f(y, y') depend on the time- function t en- 
tering into the definition of Gk, € , but for k > 2, this dependence vanishes in the 
coincidence limit y = x = y'. Similarly, Wj^La: and V Vb' [ x are indepen- 
dent of k for k > 2. 

(7) One purpose of using the point-split rcnormalization by subtraction of the 
symmetrized Hadamard parametrix % is that the latter is a locally constructed 
geometric quantity which is state independent, so that : <j> 2 : shp and (5 a b : (f> 2 : shp 
become local, covariant fields. If one applies this technique to the linear scalar field 
4> on Minkowski spacetime, one finds that : <p 2 : SH p and 3 M „ : <f> 2 : shp deviate 
for to > from the usual flat-space definitions of : (j) Q 2 : and 9 Miy : cf) 2 :, described 
above, by constants. Concretely, using the expression for Wlf*" — Sffc [k > 2), one 
can calculate that, for each Hadamard state lo of <f> , one has 

(18) lo (- <\>o 2 ■ shp{x )) = (f) 2 : (x a )) + c , m , 

(19) uj (B^ : (j) 2 : SH p(a; )) = w G (3^ : cj) 2 : {x a )) + c 2 , m V^ , 

at all points in Minkowski spacetime (details in Appendix B). Here, to is the mass 
parameter of the linear scalar field, and n^ is the Minkowski metric. The constants 
co. m and C2, TO vanish for to = 0; for to > 0, they are given by 



(20) C , m = 



(21) 02,.. 



(4tt) 2 
,4 



In ~— - 1 



TO' 



In 



(4tt) 2 

where 7 denotes the Euler-Mascheroni constant. This needs to be taken into ac- 
count in the definition of thermal equilibrium states below. 

We can now define the concept of a local thermal equilibrium state on a globally 
hyperbolic curved spacetime (M,g). Let e = (eo, ei, e 2 , e.3) be an orthonormal 
tetrad at x € M, with eo timelike and future-pointing. Then e induces an iden- 
tification of T X M with Minkowski spacetime M a , whereupon e is identified with a 
basis of M a , again with eo timelike and future-pointing in Minkowski spacetime. 
This identification is used in the following definition. 

2.3. Definition. Let lo be a state with two-point function of Hadamard form for 
the quantized linear scalar field <fi on a globally hyperbolic spacetime (M,g). Then 
let <p denote the quantized linear scalar field, with the same parameters as (f>, on 
Minkowski spacetime M a . 

(2) 

(a) We say that lo is S x -thermal at a point x G M if, with some orthonormal 
tetrad e = (e , ei, e 2 , e 3 ) at x such that e is timelike and future-pointing, there is 
a thermal equilibrium state uo^ e of <p D so that - upon identification of e with a basis 
tetrad of M Q - the equalities 

(22) w(: 4> 2 ■ SH p(aO) = <4 e (- ^ ■ shpK)) 

= Ji e {: 4> 2 : (x )) +c , m , 

(23) v a w b u;(d ab : 4? ■ shpW) = «V^ e (5V : <t>o 2 ■ shp(x„)) 

= vVuj^^ : </> 2 : (x a )) + c^w"^ 

hold for all (spacelike} vectors v,w € T X M with coordinates v^e^ — v, w v e v = w, 
for some x D e M a . (By translation-invariance of co^ e , the particular choice of x D 
is irrelevant.) 

(2) 

(b) Let N be a subset of M. We say that lo is S y N -thermal if 
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w(: <fi 2 ■ shp(^)) and to(Q a b ■ 4> 2 ■ shp(#)) are continuous in x £ N and if, for each 

(2) 

x € N, u) is S x -thermal at x. 

(2) 

The definition of S x -thermal states demands the coincidence of expectation 
values of the SHP Wick square of 4> and its balanced derivatives up to second order 
with the thermal equilibrium situation in flat spacetime. This amounts to saying 

(2) 

that Sx consists of linear combinations of the unit operator 1 and of the quadratic 
forms : <fi 2 : shp(^) and S fc : <ft 2 ■ shp(^) whose evaluations (i.e. expectation values) 
on states u> are given by (fTB|) and (fT?)) . Thus, for the linear scalar field on Minkowski 

(2) (2) 

spacetime, Si is a small subset of S x , and thus an S x -thermal state fulfills less 
constraints on its thermal properties than an -thermal state. We shall not follow 

(2) (21 

up that matter at this point. Our definition of S x -thermal states (or Sq -thermal 
states) turns out to be sufficient to derive quantum energy inequalities. 

(2) 

Given an S x -thermal state uj, we shall now use the abbreviations 

(24) iT(x) = _,(: cf 2 : SHP (x)) , 

(25) e%,{x) = -^(3 ab :0 2 : S hp(x)). 
Then we have 

2.4. Lemma. Let <f) be a linear scalar field on {M,g), with mass parameter m, and 

(2} _____ _____ 

let lo be an S x -thermal state {x € M), satisfying (|22p and (|23[) for some (3 > 
and an orthonormal tetrad e = (eo, e±, e2, 63) at x. Then the following statements 
hold. 

(a) d u (x) = jpX0,m(@) + Co, m , with 

p 2 dp 



( e Vp 2 +/3 2 ™ 2 - l)y/ p 2 + /pm 2 

(b) e u a a (x) = m 2 Xo ,m (/?) - c 2 , m 

(c) Suppose that v is a lightlike vector at x, v a v a — 0, or a timelike vector at x 
with unit proper length, v a v a — 1, and set v° — (eo) a v a . Then one has the bound 

( 26 ) « 4 )^T " v « vaCj f ~ * - v a v aC2 '' 

where 



(27) X2,m(/3) = ^ f 



p 2 J p 2 + [3 2 m 2 , 
dp 



eV p 2 +/3 2 ™ 2 _ \ 
(C(4) is the value of the (-function at A.) 

Proof. The proof is based on the fact that, with respect to coordinates induced by 
the basis tetrad e, 



where (p^)p=o,...,3 = (po,p) and po — ^\p\ 2 +m 2 in the integrals. The stated 
relations then basically result from transforming the integrals into spherical polar 
coordinates. For the upper bound in (|26[) notice that the integrand is given by 



(v°^\p\ 2 +m 2 - v ■ p) 2 / (VVH 2 + m2 - l)^|p|2+ m 2 ) 
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which is bounded above by 



H \v°- ~ V - J /(e^lEl-l). 
|p| 2 + rr 



Upon integration over p, this can be bounded by the integrand 2 (v° ) 2 |p | / (e^ 3 — 1 ) , 
using that (v ) 2 — |w| 2 = 1 in the timelike case and (v ) 2 — \v\ 2 = in the lightlikc 
case. □ 

Now we introduce a set of states whose local temperature is bounded above by 
some fixed value. 

2.5. Definition. Let /3' > 0, x € M. Then we define Lpi(x) as the set of all 

(2) 

S x -thermal states u> of the linear scalar field on [M, g) so that 



(28) ^(x)<j^ X o, m (P') + c ,r 



( 2) 

If N C M, we define Jlpi(N) as the set of all S N -thermal states of the linear scalar 
field on (M, g) so that (|28|l is fulfilled for all x G N. 

In other words, u> is in Lpi(x) if the relations ([2^)1 and (f2U)l are fulfilled for 
1/(3 <1/P'. 

Now let N be an open subset of M, and let 7 : [to, n] — > iV, r 1— > 7(t) be a 
geodesic with affine parameter t, and denote by v a — 7" the tangent vector field 
of 7. By the geodesic equation, it holds that 

(t)VV a VnJ w )( 7 (T)) = ^if(7(r)). 
Consequently, we obtain for cj 6 Lp>(N) and 77 6 C 2 ((t ,ti)), 



(29) I / r ? (T)(t; a « f> V a V b ^)(7(T))dT| = I ] l 1 "{T)d"( 1 (T))dT\ 

< \w\w 

Here, 77" is the second derivative of n. 



3. Quantum Energy Inequalities 

Let (M, 5) be a globally hyperbolic spacetime, and let <p be the classical linear 
scalar field on (M, g) with mass parameter m > and conformal coupling parameter 
^. If ip is a field configuration, i.e. a smooth solution to the field equation @, then 
the corresponding classical stress-energy tensor is a (°) co-tensor field given 
by 

= (V a v(x)){VMx)) + \9ab{x)(m 2 ^(x)-{V c ^){V cV )(x)) 

H(g ab (x)V c V c - V a V 6 - G ab {x))Lp 2 {x) , a: € M , 

where G a b = Rab — \gabR is the Einstein tensor. 

Now let <p be the quantized linear scalar field on (M, g) , corresponding to the 
choice of parameters m and £. The definition of the renormalized expectation value 
of products and derivatives for the quantized linear scalar field <f> in a state u) having 
two-point function of Hadamard form proceeds, similarly to what was done in the 
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previous chapter, by point-splitting and subtraction of the SHP (see [3^1 HI]). To 
this end, we define: 

w(: 0V„</> : sin »(x)) = V a ,W^ P Lx 

w(: 0V Q V b : SH p(a;)) = V * V* U 

w(: (V Q 0)(V 6 0) : SH P (a) - V V V U 

with k > 2, x £ M . (Note again that a and a' are identified upon taking the 
coincidence limit y = x = y' on the right hand side of each equation.) Owing to 
the symmetry of (k > 2), one can easily check that the following Leibniz rule 

is fulfilled for SHP Wick-products involving derivatives: 

w(V n (: cj> 2 : shp (a))) = 2w(: (fN a cj> : shp(#)) , 
w(V a (: 0V 6 : shp(x))) = w(: (V a 0)(V 6 0) : S hp(^)) + w(: (jN a V b cj> : sHpfa)) • 

The renormalized expectation value of stress-energy is then obtained via replacing 
the classical expressions (fi 2 {x), (V a (p(x))(X7b<p(x)), and so on, by ui(: (f> 2 : S hp {%)), 
ui(: (V a 0)(Vi0) : shp(x)), etc. Using also the Leibniz rule for SHP Wick products, 
this leads to 

^(r a s fc HP (^)) - </>V a V b : SHP (a;) + Jv a V b : </> 2 : SH p(x)) 

+ Q - A w (V Q V 6 : 4? : sw {x) - g ab {x)V c V c : <t> 2 : sw {x)) 

+ ^9ab(x)uj (: <j)V c V c 4> : S hp(i) + m 2 : (j) 2 : SH p(a;)) 
-S,G ab (x)uj(: <f) 2 : sw {x)) 

This expression, however, has the defect of a non-vanishing divergence. The way 
to cope with this problem, following Wald [35] IHj> is like this: It can be shown 
that \7 a uj(T^ b HP (x)) — VbQ{x) 1 where (apart from a free constant which can be set 
to a preferred value depending on the mass parameter m) Q is a function which is 
determined by the local geometry of (M, g); in particular, Q is independent of the 
state u). One may therefore subtract the term Q(x)g a b(x) from Lo(T^ b HP (x)) to make 
the resulting quantity have vanishing divergence. There remains an ambiguity in 
that one may still add other (°) co-tensor fields C a b which are determined by the 
local geometry of (M, g) and have vanishing divergence. We take here the same view 
as put forward in [18j , namely that the specification of C a b is a further datum of the 
underlying quantum field <fi on (M,g), in addition to the parameters m and £. An 
alternative, elegant method has been proposed by Moretti [55], which nevertheless 
we won't follow here mainly because we would like to maintain close contact to 
other works on quantum energy inequalities. This understood, we finally define the 
renormalized expectation value of the stress energy tensor in some state cu (with 
two-point function of Hadamard form) of the linear scalar field <fi on (M, g) as 

(30) u>{TT(x)) = cj(T^ p (x)) ~ Q(x)g ab (x) + C ab (x) , x e M . 

Again note that Q and C a b are state-independent and constructed locally out of 
the spacetime metric g = g a b- 

Let us next observe that, for each state u> of <j> with two-point function of 
Hadamard form, and with R — denoting the scalar curvature, 

F(x) = ^:0(V a V a +TO 2 +6R)0: SH p(x)) 

= w(: <KV Q V Q 0) : shp(z)) + (m 2 + £R{x))l>(: cj) 2 : SHP (z)) 

is a continuous function of x G M, independent of the state uj, entirely determined 
by the local geometry of (M, g) and the parameters m and £ of <f>. To see this, note 
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that 

w(: 0(V a V a + m 2 + ZR)4> : SHP (z)) - (V a V a , + m 2 + £R)W^U . 

On the other hand, W^ P (j/, 2/') is the integral kernel of 1TJ| P = ~%(k>2), 
and since #^(/, (V b \7 b + m 2 + £R)h) = 0, it follows that (V a V a , +m 2 + £i?)W^ p 
is independent of lu as the w-dependent part of W^F is annihilated by the wave- 
operator (V a 'V a <+m 2 + £i?). Inconsequence, F(x) = (V a/ V a / +m 2 + £R)W^f [ x 
is state-independent, continuous in x, and actually it is determined by the local 
geometry of (M, g) since so is (by the Hadamard recursion relations) . 

Using the Leibniz rule, we can now rewrite the expression for uj(T^ n ) as follows: 

(31) to(T^ n (x)) = w(- : 0V a V b : SHP {x) + iv a V h : 4> 2 : SHP (a;)) 

+(- -C)w(V a V b : <j> 2 : BB p(ar)) 

+ (4£ - 1) M- : 0V c V c : SHP (x) + ^V C V C : 4> 2 : SH p(aO)) . 9a6 

+ (Y(l - 40 (™ 2 + £#) - \iR\ 9ab - £G o6 ) w(: 2 : SHP (z)) 

+ ((^ - 5) W - + C afc (a;) 

By (|15[) and once more the Leibniz rule, e ab can also be expressed as 
fab = w(- : ^VaVb'/' : shp(z) + ^V a V fc : 2 : sHpfa)) 

Thus, if u> is an S x -thermal state of <p, we obtain 

(32) w(KF(x)) = e&(x) + (4£ - l)g ab (x)e» c (x) 

+ (\ - OVaVbiT (a;) + (gafc(x)V(x) - ^(a;))^ (z) 
((4£ - l/2)F(a:) - Q{x))g ab {x) + C ab {x) , 
where we use the abbreviation 

(33) ip{x) = (1 - 4£)(m 2 + £R(x)) . 

With this expression, we are now in the position to derive bounds on »°w w(T™) 
for lightlike or timelike vectors v. We will treat lower bounds first. 

3.1. Theorem. Let <fi be the quantized linear scalar field on (M, g), with parameters 
m,£ and C a b, o,nd let U be a state of '</> having two-point function of Hadamard form. 

(a) Suppose that £ = 1/4, and let v be a lightlike vector at x £ M , or a timelike 
vector at x with v a v a = 1. If to is in Lpi(x), ft 1 > 0, then 

(34) v a v b w(T% n (x))> q{x,v;0) 
where 



q{x,v;(3') = -- \v a v b R ab (x)\ 



-^XO.mC/?') + CO,- 

p 



^ - Q(z) - ic 2 , m ) Wa i; a + </VC a6 

(b) Let £ 6e arbitrary, let N C M, and Zet 7 : [to,ti] — > N be an affinely 
parametrized lightlike geodesic defined on a finite interval, with tangent 
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vector field v a — 7 . Suppose that r\ is in Cq((tq,ti)) with 77 > 0. If 
u> E £jp>(N), there holds the bound 



(35) 



Here, writing 



R [7] = max | 7 a (r)7 6 (r)i? at ( 7 (r))| , 
anrf defining C[ 7 ] analogously, the bounding constant is given by 



l£|-R[7] 

r — -ei - 



-^2X0,m(/3') + C , m 



•c. 



[7] 



(c) Let £ be arbitrary, let N C M, and let 7 : [to,ti] — > N be an affinely 
parametrized timelike geodesic with tangent vector field v a = j a , so that 
v a v a = 1. Assume that n is in Cq{{tq,ti)) with r\ > 0. If w E Lp/(N), 
there holds the bound 

(36) J V (T)v a v b uj(T:r(l(r)))dT > Ql (7,^/3') , 

where, using the notation ipLyi = max T£ [ ToTl ] |V'(7( r ))l; an d defining -Fr 7 i 
and Qw\ similarly, the bounding constant is given by 



gi(7,r?;/3') = — I ^ — CI 



-2X0, m(/3') + C ,m 



l^'IlL! 



2XO,m(/3') + Cq,. 



1*7 Li 



-(V^ + lel^]) 
- ^ - i|F [7] + |4C - l||c 2 , m | + Q h] + C h] + 



C2,-, 



l»7 i» 



Proof. The proof of the statement consists just of inserting the estimates of 
Lemma 12.41 and discarding manifestly positive terms, in combination with estimate 
([29| for the average of the second derivatives of along the geodesic. The term 
involving second derivatives oft?" doesn't occur for £ = 1/4, which makes it possible 
to give a pointwise lower bound in this case. □ 



The central assertion of Theorem 13. II is that the lower bound of the energy density 
averaged along a causal geodesic depends only on the temperatures an LTE state 
attains on the geodesic, and is otherwise state- independent. The bound worsenes 
(shifts towards the left on the real axis) as the temperature increases, i.e. with 
increasing 1//3'. It should also be noted that the bounds are local covariant. 

For upper bounds on the averaged energy density of LTE states, an additional 
state-dependence shows up: The bounds depend also on the tetrad e appearing in 
the condition of Sx -thermality, Def. 12.31 In this sense, the lower bounds on the 
averaged energy densities of LTE states are stronger than the upper bounds. This 
is similar to what holds for averages of energy densities for arbitrary Hadamard 
states of the linear scalar field [2] . 

Let x E M, and let e = (eo, . . . , 63) be an orthonormal tetrad at x with eo timelike 
and future-pointing. We define Zpi(x,e) as the set of all states u> in Zp/(x) where 
the S x 2 ^ -thermality conditions (|22[) and (|2"3"|) hold with respect to the given tetrad. 
Similarly, let N be a subset of M, and let N 3 x ^ e(x) = (eo(x), . . . , e^ix)) be a 
C° field of orthonormal tetrads over N, with cq(x) timelike and future-pointing for 
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all x. Then we define &pi(N, e) as the set of all states u) in Lpi(N) such that, for 

(2) _____ _____ 

each x £ N, uj satisfies the S x -thermality conditions (f2"2"| and (l2l?l) with respect 



-^2Xo,m(/3') +c . 



to e = e(ic). With these conventions, we obtain the following upper bounds on 
(averaged) energy densities. 

3.2. Theorem. Let <f> be the quantized linear scalar field on (M, g), with parameters 
m, £ andCab, and letuj be a state of<f> having two-point function of Hadamard form. 

(a) Suppose that £ = 1/4, let v be a lightlike vector at x € M, or a timelike 
vector at x with v a v a — 1, If u> is in Lp>{x,e), f3' > 0, then 

(37) p( V ;J?',e)>// W (CW) 
where 

p(v,x;p',e) =C(4)^J + \ \v a v b R ab 

+ q{x,v;0) 
with v° = v a {eo) a - 

(b) Let £ be arbitrary, N C M, and let 7 : [to,ti] — > N be an affinely 
parametrized lightlike geodesic defined on a finite interval, with tangent vec- 
tor field v a — 7 . Suppose that r\ is in Cq{{tq,t\)) with 77 > 0. If ui is in 
Lp>(N, e), then 

(38) P0 ( %m p',e)> J V (r)v a v b uj(T r J n h(r)))dr 
where 

Pod, m 0, e) = ^L( v ^\\rj\\ L1 + | 50 ( 7 , /?'; ^ 
with u° 7] = max Te[r0!ri] 7o(T)eg(7(r)). 

4. Averaged Null Energy Condition (ANEC) 



In this section we derive the averaged null energy condition (ANEC) for S 



(2) 
N 

thermal states of the quantized linear scalar field on a globally hyperbolic space- 
time (M, g). 

The ANEC on a state ui of <f> demands that 

(39) liminf / + v a v b Lo(T r a ^ n { 1 { T ))) dr > 

for all complete geodesies 7 in M with affine parameter r and tangent v a = 7 . 
If this condition holds, and if (M,g) together with tj> and u> are a solution to the 
semiclassical Einstein equation in the form 

(40) G ab (x) = 8%u(T% n ( x )) , x e M , 
then this implies that 

(41) liminf / + v a v b G ab (~f(T))dT > 

for all complete geodesies 7. (We address the issue for the semiclassical Einstein 
equations with an additional contribution by a classical stress-energy tensor below.) 
It has been shown that this weaker form of the usual pointwise null energy condition, 
which demands that £ a £ b G a b(x) > for all lightlike vectors l a at each x € M, is 
still sufficient to reach the same conclusions with respect to singularity theorems as 
obtained from the pointwise null energy condition, i.e. that congruences of geodesies 
will focus with expansion diverging to —00 at finite affine geodesic parameter |25] . 



18 



The validity of (|39p is therefore of importance for the properties of the spacetime 
structure of solutions to the semiclassical Einstein equations. 

It has been argued in [42] that condition (|39|) may be replaced by the following 
condition: 

/oo 
f]x{T)v a v b uj(T^ n (-f{T))) dr > 
-oo 

for any i] <E C%(R), rj > 0, with r)(0) > and T) X (r) = r?(Ar) for A > 0. More 
precisely, in [42] it has been shown that |42|) and (f40|) imply that the expansion 
of a congruence of lightlike geodesies around 7 becomes singular along 7 (in the 
sense of diverging to — 00 at a finite value of the afhne parameter) unless it vanishes 
identically on 7. (In [12] this argument is given for half-line geodesies, but it carries 
over to the case at hand as will be shown in our Appendix A.) 

(2") 

Now let to £ S N , and let 7 be a complete lightlike geodesic in N C M with 



affine parameter r and tangent v a — j a . Then, from (1321) . 

(43) «Va,(T Q 7;") = v a v b £ab + (j ~ £) v a v b V a V b #" - £v a v b G ab $ u + v a v b C ab 

holds along 7. Therefore, positivity properties of the (integrated) energy density 
v a v b uj(T^ n ) depend also on the behaviour of G ab and C ab . The sign of the term 
involving G ab is not known. To circumvent this difficulty, we assume that the un- 
derlying spacetime (M, g) together with (f> and u> are solutions to the semiclassical 
Einstein equations (|40[) . since it is this situation in which the ANEC is applied to 
deduce flTj) and the ensueing statements about focussing of lightlike geodesies. Sup- 
posing that [M , g) together with <j> and u> are solutions to the semiclassical Einstein 
equations, and also that u> is an -thermal state, we obtain upon combination of 
(|3D| and |03]) the equation 

(44) v a v b [G ab (l + 8tt^) - 8irC ab ] = 8irv a v b (e w ab + ( j - A V a V b ^ 

on N. In order to draw further conclusions, one must specify C ab . We recall 
that C ab is a datum of the linear quantum field (f>, a priori only restricted by the 
requirement that T^§ n be a local covariant quantum field and divergence-free, thus 
C ab should be locally constructed from the spacetime metric. Following Wald [4Tj . 
one can make the assumption that C ab have canonical dimension, which leads to 
the form 

(45) C ab = Ag ab + BG ab + T-^S^g) + D^- b S 2 {g) 

where S±(g) — J M R 2 dvo\ g , 82(g) — J M R ab R ab dvol g , and S/Sg ab means functional 
differentiation with respect to the metric, with constants A, £?, T, D as remaining 
renormalization ambiguity for the quantum field <f> (see |41] for additional discus- 
sion). For the rest of our discussion, we will simplify matters by assuming T,D = 0. 

Making this assumption, so that ([43)) holds, and observing that hence, v a v b C ab = 
Bv a v b G ab for all lightlike vectors v a , (|4"4")) assumes on N the form 

(46) v a v b G ab (l + 87r(fiJ w - B)) = ^v a v b (e u ab VaV 6 #* 

The constant B is still free, and one may now try to choose B in such a way that 

(|4"6")l entails the ANEC for all lightlike geodesies in N C M and an as large as 

(2) 

possible class of S N -thermal states u>. We will show that this is possible with 
different conditions on B for the cases ^ = l/4,0<^< 1/4, £ = 0. 



If) 



4.1. Theorem. Let (M,g) be a globally hyperbolic spacetime, let (f> be the quantized 
linear scalar field on (M, g), with parameters m, £, C a b, where C ab — Ag a b + BR a b, 
with real constants A,B. 



(2) 

Suppose further that ui is a quasifree Hadamard state for cf>, that u> G S N for N 



C 



M, and that (M, g) together with ip and uj provides a solution to the semiclassical 
Einstein equation (|40[l . 

Let j be a complete lightlike geodesic in N with affine parameter r and tangent 
v a = 7°, and let r) G Cg(K), r\ > 0. Then 



(47) Urn/ n{XrMT:r{l{r)))dT>Q 

A ^ u J -oo 

holds if any of the following groups of conditions is assumed: 

1. ) £ = 1/4, B < 1 + 27rco, m . In this case one even has 

v a v b u;(T% n (x)) > 
pointwise for all x G M and all lightlike vectors v a at x. 

2. ) < i < 1/4, B < £c , m + l/(8rr), 

(48) Aln(i? u, (7(r/A))) -► as A^O for almost all r, 

(49) ^ A| ln(iT (7(t/A)))| dr < fc < oo for small A and s < r G 

3. ) £ = 0, S < 1/87T, 

(50) AiT(7(t/A)) -> as A^O for almost all r, 

(51) / X^(j(r/X))dT < K < oo for small X and all s < r G K . 



Remark (a) If, instead of ppf . the semiclassical Einstein equations are 
assumed to hold in the form 

G ab (x) = te(T$>"(x)+u>(T%»(x))) 

with a stress-energy tensor T£ l b ass for classical, macroscopic matter distribution, 
and if it is assumed that this stress-energy tensor fulfills the pointwise null energy 
condition e a £ b T^ l b ass (x) > for all lightlike vectors i a at each point x G M, then the 
statements of the theorem remain valid with Tf b ass + uj{T^ b n ) in place of uj(T^ b n ). 

(b) Conditions (|4"8"|) and ([41?]) say, roughly speaking, that $ w (7(t)) shouldn't grow 
faster than e |r|<1 " e) for |r| -> 00, while ([50]) and ([51]) say that ^(7(7-)) shouldn't 
grow faster than |r| 1_e as |r| — > 00. Now since $"(7 (r)) = (/3(7(r))~ 2 xo,m(/3(7(T)))4 
co, m and since 

1 Z" 00 

Xo,m(/3)^^2 j f p dp for /3^0, 

this means that the growth of the temperature 1/ (3{^{t)) at 7(7") appearing in Def. 
12.31 of S^L-thermality u should not exceed e' T ' a e> 2 and |r|( 1_e ^ 2 as |r| — > 00, 
respectively. 

Proo/ 0/ TTim. 1.) If £ = 1/4, then (jUJ) assumes the form 

(52) uVG a6 (l + 8tt(#74 - B)) = 8irv a v b e% b . 

If B < 1 + 2irco. m , then the factor 1 + 87r($ w /4 — B) is strictly positive, as is the 
right hand side of ([52"]) . This equality holds pointwise at all x E M and for all 
lightlike vectors v a , thus proving, in combination with the assumed property (|40[) , 
the statement of the theorem. 
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2. ) For < £ < 1/4, B = £c , m + 1/(8tt) - £c, where c > 0, gSJ) takes the form 
(53) «VG afc (87r£(7r - c , m + c)) = 8irv a v b e% b + 8tt(1/4 - £K« b V a V b ?r . 
Observing that w a t |b V Q VfcCo l m = 0, the last equation is turned into 

at r v a v b ^ ab , (l/4-Q^ 6 V a V 6 (^-co, m ) 

( ' VVG *»~ - Co>m + C ) + £(#o - c ^ m + C ) 

where it was used that d u — co, m + c > 0. The first term on the right hand side of 
(|54|) is positive. Upon integration against a non-negative Gq weighting function -q 
along the geodesic 7 we obtain, using the abbreviation 

u(r)=^( 7 (T))-co, m , 

the inequality 

^(r)K W b G afc )( 7 (r))dr > ^- £ [ v (r)^f^-dr. 

£ J u{t)+c 

By partial integration, 

u(t)+c J \u(t) + cJ J 

Thus, since the first integral on the right hand side is non-negative, (1/4 — £)/£ > 
for the £ considered and using the monotonicity of the logarithm together with 
c > 0, 

ry(Ar)K« b G ah )( 7 (r))dr > ±l±zA J Aln(u(r / '\))ff' *(r) dr , 

and owing to assumptions (|4"5|) and |4^|) . the expression on the right hand side 
converges to as A — > 0. Equation (|4T|l is then again implied by the assumed 
property ([40]) . 

3. ) If £ = 0, equation (|46p turns into 

(55) w°w 6 G a6 (l - 8^5) - Ww 6 ^ + ^VV a V 6 <r , 

and by the condition on B, the factor 1 — 8irB is strictly positive. Observing 
again positivity of 8wv a v b £ J b , upon integration against a non-negative Cq weighting 
function 77 along 7 one obtains 

r,{\r)v a v b GMr))dT > | A U (t/AV'(t) dr 

and the right hand side converges to as A — » by assumptions (j5U|) and ([STj) . 
Again (j4"T|) is deduced from the assumed validity of (JlQ"]) . □ 

5. Generalized Local Thermal Equilibrium States 

_ (2) 
The notion of LTE states in [9J, and the related definition of S x -thermal states, 

is actually more general than the definition given in Sec. 2. In [5j the possibility 

was considered that an LTE state w coincides on S^-observables not necessarily 

with a thermal equilibrium state at sharp temperature in a certain Lorentz frame, 

but with a mixture of such states. 

In our setting, where we work with the linear scalar field, this corresponds to a 

modification of Def. 2.3 as follows. As a consequence of eqn.(fl"2"|). ui^ e , the quasifree 

thermal equilibrium state with respect to the Minkowski tetrad e = (eo, ei, e2, ea) 

at inverse temperature (3, depends only on (3 = /3eo. This quantity completely 

parametrizes w^ e , so we write u)@ in place of ui^ e . The vectors (3 take values in V + , 

the set of future-directed timelike vectors in Minkowski spacetime. 
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Let (M,g) be a globally hyperbolic spacetime, let V x C T X M be the set of 
future-directed timelike vectors at x £ M, and let p x be a Borel measure on V x 
supported on a compact subset B x C V^ + , with J B dp x ((3) — 1. Then we say that a 

( 2) 

Hadamard state w of the linear scalar field <fi on (M, g) is a generalized S x -thermal 
state if 

w(: 2 : shp(s)) = / ( : 0« 2 : {x ))dp x (f3) + c Q . m , 

v a w h Lo(d ab : 2 : shp(x)) = vV f «?(3^ : <^ 2 : (s o ))dp09) + c 3 , m v"ti; , V 

holds for all (spacelike) vectors v, w £ T X M for some x £ M Q . Making further the 
assumption that F h-> / M J* B F(x, f3)dp x {f3)dvo\{x), F £ C (TM, C) is a distribu- 
tion (on the manifold TM), such that x \— > J" F(x, f3)dp x ((3) is C 2 , one can define 

(2") (2) 

generalized -thermal states in analogy to the definition of -thermal states 
in Sec. 2. 

With these conventions and assumptions, the results of Thms. 3.1, 3.2 and 4.1 

(2) 

extend to generalized Sjf -thermal states, under identical assumptions, except that 
the bounds have to be corrected for the p^-integrations. It should be obvious how 
this is to be done. 



6. Discussion and Outlook 
We have generalized the concept of local thermal equilibrium states of [S], or 

(2) 

rather, the concept of Sx -thermal states, to the quantized linear scalar field mod- 
els on generic globally hyperbolic spacetimes, and have shown that one can derive 
certain quantum energy inequalities for such states. The lower bounds appear- 
ing in the quantum energy inequalities of local thermal equilibrium states depend 
only on the local temperature of the states, i.e. thermal function corresponding 
to the expectation value of the Wick-square in local thermal equilibrium states. 
The upper bounds, instead, depend also on the local frames with respect to which 
S x -thermality is defined. In this sense, the lower bounds are stronger (have less 
dependence on the states) than the upper quantum energy inequality bounds. This 
is a feature also found for quantum energy inequalities of general Hadamard states 
of the linear scalar fields, and has led to the proposal to consider 'relative quan- 
tum energy inequalities' as a more general variant of quantum energy inequalities 
which has the potential to be valid also in interacting quantum field theories [TH] , 
Moreover, the quantum energy inequalities for thermal equilibrium states are local 
covariant. 

A major purpose of quantum energy inequalities, especially in local covariant 
form, is to provide information about the structure of spacetime geometries ap- 
pearing as solutions to the semiclassical Einstein equations. Quite generally, they 
serve as stability conditions on quantum matter, and ensure that correspondingly 
the (semiclassical) gravitational interaction is attractive, at least when averaged 
over sufficiently extended spacetime regions. The averaged null energy condition 
which we proved for certain values of the curvature coupling £ and certain values of 
the renormalization constants is of a similar nature. One may also take the require- 
ment that the ANEC should be fulfilled for suitable thermal equilibrium states as 
a constraining condition on the largely free choice of renormalization constants for 
the stress-energy tensor. Certainly a demand in this spirit leads to further rela- 
tions between the renormalization constants, the parameters fixing the field model, 
and possibly geometrical quantities, and for this reason it is attractive to further 
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study quantum energy inequalities and ANEC in the context of solutions to the 
semiclassical Einstein equations. 

One important issue we haven't addressed at all so far is the existence of local 

(2) 

thermal equilibrium states, or at least S\r states for subsets N in spacetime. We 
have simply assumed that there are such local thermal equilibrium states to which 
our results apply. The question if there are local thermal equilibrium states in 
generic spacetimes is an interesting and difficult problem, for which we can't offer, 
as yet, any route to its solution. However, the existence of LTE states for the 
massless and the massive Klein-Gordon fields on (parts of) Minkowski spacetime 
has been established, with an interesting relation to situations resembling a big 
bang scenario [H [27J . 

The question if local thermal equilibrium states exist is a first step towards the 
question how generic they are. One is inclined to think that within certain time- 
and energy scales, local thermal equilibrium states should be the archetypical phys- 
ical states in the sense that, if one is asked to randomly pick a state in the physical 
state space (of a quantum field theory), then the result would be a local thermal 
equilibrium state with overwhelming likelihood. At least this is expected for in- 
teracting quantum field theories since interaction tends to equilibrate subsystems 
(or degrees of freedom) of a large system. If this turned out to be true, and if the 
temperature distribution of such states turned out to allow ANEC results similar 
to those of Thm. 4.1, then one would be led to conclude that (under general addi- 
tional assumptions) the occurrence of singularities in solutions to the semiclassical 
Einstein equations is a generic feature. It would be of utmost interest to investigate 
this circle of questions further particularly in scenarios of early cosmology. 

Appendix 

A We will present a result on real- valued solutions 9(t) of the differential equation 



It follows from the Picard-Lindelof Theorem that there is an open interval (a, b) 
containing 0, which may be finite, semi-finite or infinite (i.e. coinciding with R), 
such that this interval is the domain of the unique, inextensible C 1 solution 9 of 
([56]) satisfying the initial condition. In this case, we call 9 the maximal solution of 
(|56|) defined by the initial condition, and refer to (a, b) as the maximal domain. 

The following statement is a variation on a similar result in [32], and it uses 
a very similar argument, the main difference being that the assumption ([55]) here 
is slightly different from that in [12], where the integral is taken over a semi-axis. 
Note also that our parameter A corresponds to 1/ A in the notation of [32] ■ 

Theorem A.l Suppose that f £ C 1 (M, R) has the property 



for the function T](t) = (1 - t 2 ) 4 for \t\ < 1, r/(t) = for \t\>l. 

Then either the maximal domain of 9 coincides with all of the real axis and 
6 {t) — for all t G R, or the maximal domain (a, b) of 9 is a finite or semi- 
finite interval. In this case, 9(t) — > ^Foo for t approaching the finite boundary at 
the right/left side of the maximal domain (in the finite case this holds with the 
respective sign for both boundaries). In particular, this is the case if 8 (to) ^ for 




(58) 
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some to in the maximal domain of 9. 

Proof. Consider the auxiliary differential equation 

(59) u"{t) + —u(t) = 

A 1 

For the initial values u(0) = 1, tt'(O) = Oq and the given / this linear differen- 
tial equation has by the Picard-Lindelof Theorem a unique, global solution u G 
C 2 (K, M). Furthermore, this solution is nonzero in some neighbourhood of 0. For 
points from this neighbourhood, one can then rewrite (|59[) as 



d (u'{t)\ , (u'(t)\ 2 f(t) 



dt V u(t) ) V U M / U M 

which implies that §(t) = fulfills equation (I56|) . Furthermore, 6 also satisfies 

the initial condition ([57)) and by the uniqueness part in the Picard Lindelof Theorem 
it therefore agrees with 9. This however implies that the only way in which 9 can 
fail to be C 1 at a boundary point c = a or c — b of a semi-finite interval is a zero of 
u at lie. At this zero v! has to differ from zero, otherwise u as a C 2 -solution to (l5§|) 
with initial conditions u(pc) — 0,u'(pc) = lim^^pc u'(x) = would be identically 
zero in contradiction to the initial values for u at 0. By continuity, u' is therefore 
nonzero in a neighbourhood of /ic, and by (|56|) . 9(t) = approaches the value 

—oo for t — > c, t < c (right boundary point) or the value +00 for t — > c,t > c (left 
boundary point). For proving that 9 diverges at the boundary (boundaries) of a 
semi-finite interval it is therefore sufficient to show that 9 cannot be continued as 
a C 1 function beyond this boundary. 

With the definition of r\ as above, and provided that the maximal domain of 9 
coincides with all of R, one has for < A < 1, 

/oo 
8(t)Xr)'(Xt)dt 
-OO 

,1/A 

= -8A / 9{t){Xt){l - {Xtffdt 
J-i/\ 

> -8A / \9{t)\ (I -(Xtffdt 
J-i/x 

owing to the fact that both \Xt\ and |(1 — (Xtf)\ are bounded by 1 on the domain 
of integration. Combining this with ([55]) and |55|) leads to 



,1/A ,1/A 

(60) limsup-8A/ \9(t)\(l - [Xtff dt + fi / 9(tf(l - (At) 2 ) 4 dt < . 

A^0 J-l/X J-l/X 

Using also the Cauchy-Schwarz inequality 

l/A / ,1/A 

\9{t)\{\ - {Xtff dt<\ 9{tf{\ - {Xtff dt 

l/A V/-A 

the estimate (|60l) can be replaced by 
(61) 

o / ,1/A \ 1 I 2 ,1/A 

limsup --V2A / 9{tf {I -{Xtffdt) + 9(tf {I - (Xtff dt < , 

A^0 A 1 \J-1/X J J-l/X 

which shows that 9(tf dt — upon using Levi's theorem. Since 9 is C 1 , this 
implies that 9(t) = for all t. 
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We have therefore shown that the assumption of 9 being C 1 on all of R implies 
9{t) = for all t € R; if on the other hand 9 is C 1 only on a maximal finite or 
semi-finite interval, then by the statement in the first paragraph of the proof, it 
will diverge at the finite boundaries of this interval in the indicated way. 



B Here we will calculate the constants Co,m> C2,m that arise when defining the Wick- 
square and the second balanced derivative on Minkowski spacetime using the co- 
variant point-split renormalization. A similar calculation can also be found in the 
Appendix B of [2H], the different conventions adapted here however lead to small 
changes in some of the formulas appearing. 

The Hadamard recursion-relations satisfied by the functions Uj in (jTDJ) read with 
our sign-conventions: 

-2(W)V„E7b - (4 + V K V K a)C/ = (V K V K + m 2 + £R)U 



(V K V t 



-2{V K a)V K U j+1 + (4j - V K V K a)U J+1 = 



For Minkowski spacetime, U is identically one, V K V K cr 
tions of the resulting recursion relations 



i + i 

-8 and the unique solu- 



4(x-x') K V K U + 4U Q = m 2 



4(z - x') K V K U j+1 + 4(2 + j)U j+1 = 



(m 2 + V K V K )U 3 
J+l 



that remain bounded for x 
characteristics) as 



x' are easily calculated (e.g. using the method of 



1 



2 \ J+ 1 



J iKi + i)! 

For non-lightlike x — x' where Gk,e is a regular distribution (the corresponding 
function being obtained as the pointwise limit e — > 0) we have with the abbreviation 
{x-x'f :=r lab (x-x') a (x-x') b -- 



1 



\x - x'f 



i 



— m 2 {x — x'Y 



(it will be seen in the course of the calculation, that Gi.o is actually sufficient to 
calculate the second balanced derivative, one does not need C?2,o)- The two-point 
function Wg of the Minkowski vacuum state u> vac for spacelike (x — x') is given 
by ® 



Wf*°(x,x') 



m K X (mV '-{x-x'f) 
4^2 



^-(x-x'y 

and using the asymptotic expansion of the modified Bessel function K\ for small 
arguments, the terms up to the order (x — x 1 ) 2 of the two-point function are given 
by 



1 



-(x — x') 2 

,2 



2 



-m 2 (x — x') 2 



-i / /\2 

1-— {x-x'f 



m 



(2 7 -l) + (2 7 -5/2)- 



-m 2 {x — x'f 



(here and in the following, x — x' is now assumed to be spacelike). The difference 
Wg" {x + C; x — () — Gifi{x + C, x — () to the order required for the calculation of 
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4> 2 : shp(x)) and w°°(^^ v : 4> 2 : shpOe)) is then 



Wf & \x + C, x - C) - Gi, (z + C, a: - C) 



With this expression one calculates 

<\> 2 : shp(x)) 



2 r 

m 



(4tt) 



In ( — ]+2 7 -l 



In ^- + 2 7 - 5/2 



■m ( 



2^2 



(4tt) 2 

4 

w°°(9^ : 2 : shp(x)) 







H 4 ) 


- 1 







4 7 ' 1 2 



Co.. 



In 



2 ^m 2 \ 5 
" 2 



(47T) 

and from this one reads of the equations (|T5)) and fT§]) 
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